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Abstract It is known that the torsional rigidity for a punctured ball, with the puncture having the 
shape of a ball, is minimum when the balls are concentric and the first eigenvalue for the Dirichlet 
Laplacian for such domains is also a maximum in this case. These results have been obtained by 
Ashbaugh and Chatelain (private communication), Harrell et. al. Ifl2l . Kesavan [ 13 ] and Ramm and 
Shivakumar [ 18 1. In this paper we extend these results to the case of /j-Laplacian for 1 < p < 00. For 
proving these results, we follow the same line of ideas as in the aforementioned articles, namely, 
study the sign of the shape derivative using the moving plane method and comparison principles. In 
the process, we obtain some interesting new side results such as the Hadamard perturbation formula 
for the torsional rigidity functional for the Dirichlet /^-Laplacian, the existence and uniqueness result 
for a nonlinear pde and some extensions of known comparison results for nonlinear pdes. 
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1 Introduction 



The /5-Laplacian A„ is the non-linear operator defined as A p f — div(| V/| /J_2 V/). Let B\ be an open 
ball in R. Let Bq be another open ball whose closure is contained in B\, and is free to move inside 
B\ . Let Q =Bi\Bq- We consider the following domain optimization problems: 

i. GivenyeW ''(X2),the unique solution of the equation 



—Apii =1 in Q, 
m = on d£2, 

we are interested in minimizing the p-torsional rigidity 



E(Q):= / \Vy\ p dx = / ydx 
Jn Jn 



(1.1) 



(1.2) 



with respect to the position of the hole Bq. 
ii. Given the eigenvalue problem 



-A p u = X\u\ p 2 u inQ, 
w = on dQ. 



whose principal eigenvalue is 

Ai (Q) := inf < 



||V< 



lp (a) 



\lp(Q) 



(1-3) 



(1.4) 



we are interested in maximizing Ai (£2) with respect to the position of Bq. 

The following results were obtained, in the linear case, i.e., for p = 2, by Ashbaugh and Chatelain 
(private communication), Harrell et. al. 02), Kesavan 11131 . Ramm and Shivakumar 0181 : the tor- 
sional rigidity is minimum if and only if Bq and B\ are concentric. Also, the first eigenvalue Ai of 
problem ( 11.4b attains its maximum if and only if the balls are concentric. 

The analogues of these results for manifolds were obtained in Anisa and Aithal [2] in the setting 
of space-forms (complete simply connected Riemannian manifolds of constant sectional curvature) 
and in Anisa and Vemuri (On two functionals connected to the Laplacian in a class of doubly con- 
nected domains in rank-one symmetric spaces of non-compact type, preprint) in the setting of rank- 
one symmetric spaces of non-compact type. We extend these results, in a different direction, to the 
non-linear setting. Our main results are Theorem l6.1l and Theorem l6.2l 

The proofs in [ 1 3 , 1 8 1 rely on shape differentiation ll20ll . the moving plane method iHl fTTI and 
various maximum principles. In the non-linear case, carrying out this program involves several tech- 
nical difficulties. We develop the shape calculus for the torsional rigidity function for p-Laplacian. A 
formula for the Hadamard perturbation of the first Dirichlet eigenvalue for the p-Laplacian is given. 
This, however, is not new and may also be seen in the works of Garcia Melian and Sabina de Lis 
iflOl . Lamberti lTT4l and Ly lfl6l . For the Steklov eigenvalue this is done in Del Pezzo and Fernandez 
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Bonder |9|. Subsequently, we analyze the sign of the shape derivative. We do this by proving a suit- 
able strong comparison result. In the case of the eigenvalue problem, before this, we also need to 
prove a general weak comparison principle for the p-Laplacian with non-vanishing boundary condi- 
tion (cf. Theorem 13. U . This result is new and can be of independent interest in itself. An existence 
and uniqueness result for a nonlinear pde is required for applying this comparison principle and this 
result is also proved (cf. Proposition l4.lt . 

The Section [2] establishes notations, contains some definitions and technical preliminaries. In 
Section [3] we recall some existing weak and strong comparison principles for the p-Laplacian and 
prove an extension of a weak comparison principle. In Section|4] we prove the existence and unique- 
ness of non-negative solution for a nonlinear pde needed for an application of the comparison prin- 
ciple. In Section [5] following [20] we obtain the Hadamard perturbation formula for the torsional 
rigidity functional (11.21 i and for the first eigenvalue of the Dirichlet /5-Laplacian ( 11.41 i. Finally, in 
Section|6]we prove the main results by analyzing the sign of the shape derivatives. 



2 Prelminaries 

In this section we introduce some definitions and recall some results which will be used later on. 

Shape Derivative: Given a functional / which depends on the domain £2 (usually, a smooth 
open set in Mr) and given, a variation of the domain £2 by a fairly smooth perturbative vector field 
V which has its support in a neighborhood of d£2, the infinitesimal variation of J in the direction V 
is defined as 

mv)=*n j w- j w (21) 

f->0 t 

where £2, is the diffeomorphic image <£>, (£2 ) of £2 under the smooth perturbation of identity <P t (x) = 

(I + tV)(x). 

The shape derivative is a tool widely used in problems of optimization with respect to the domain 
as it permits to understand the variations of shape functionals (cf. Simon [ 19], [20|). 

We define : B(t) := (D<P t )-\ y(f) := |det£><£, and A(t) := y{t)B(t)B{t)* where B(t)* shall 
denote the transpose of B(t). It will be convenient to denote 7(f), B(t), B(t)* and A(t) respectively, 
by Yt, B t , B* and A t . We observe that 

D& t =I + tDV (D&t)* =I+t (DV)* (2.2) 

and so, B t , B*, A t , Yt and F t are analytic functions of t near t = 0. We record that 

Y(0) = divV (2.3) 
(B,*)'(0) = —(DV)*. (2.4) 

So, for small t, we have 



Y w 7(0)+f7 , (0) = 1+fdivV 



(2.5) 
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Also, for t sufficiently small say \t \ < ?o, there exists a constant C > such that 

|(Z)<£ r )*£|<C|£|forall£eIR". (2.6) 
Consequently, by substituting B*T] for rj arbitrary in M", we have 

|S f *T7 1 > C 1 1?7 1 for all 77 e M" . (2.7) 

Pucci-Serrin IDENTITY: We shall find it very useful to employ the extended version of the 
Pucci-Serrin identity proved by Degiovanni et. al. [8| which gives the following identity for the 
p-Laplacian. Assume that u € C 1 (£2) is a solution of the equation 

— A p u = f in £2, 
u = on d£2. 

Then, for all V G C 1 (£2) the following identity holds 

V-ndS= f diW^^-dx- [ ({DV)*Vu,\Vu\ p - 2 Vu) dx+ ( V -Vuf dx 
Jn p Jn Jn 

(2.9) 

The Pucci-Serrin identity may be obtained by using V ■ Vu as a test function in ( 12.81 ) and after several 
integration by parts whenever u G C l (£2) DC 2 (£2). However, by standard regularity results for so- 
lutions of the p-Laplacian equation, they are known to belong to only C 1,a (£2) (cf. Tolksdorff [22]) 
as the coefficients \Vu\ p ~ 2 degenerates near the critical points of u. This formula can be justified 
by regularizing the coefficient first and then passing to the limit cf. [8 1 (see also Garcfa-Melian and 
Sabina de Lis [ 1 1 and the work of Del Pezzo and Fernandez Bonder [9] for such arguments). 

A POSITIVE DEFINITE MATRIX: Define a strictly convex function r : ->• K by r(x) — ^j-. Let 
A = Dr. Then A = (Ai,A 2 , . . . ,A N ) :R N ^R N and is given by 

A(x) = \x\ p - 2 x. (2.10) 

N 

corresponds to the symmetric matrix 

!j=l 

|x| f,_2 Id+ (p — 2)\x\ p ~ 4 x<3x, which is the Hessian of the convex function F. It can be seen that (p — 
l)\x\ p - 2 and \x\ p - 2 are eigenvalues of si with multiplicity one and (n — 1) respectively. Therefore, 
for any | e W, we have 

<si%,% >>min{l,/7-l}|x| p " 2 |^| 2 . (2.11) 




(/>-!) 





du 


L 


dv 



dA 

Clearly, A e ^°°(R N \ {0}). The matrix si := -^-(x) 
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3 Comparison Theorems for the p-Laplacian 

Let £2 C M. N be a bounded domain with smooth boundary. Let B : £2 x ffi — > W. be continuous function 

and m i — > fi(x,u) is locally Lipschitz on R \ {0} uniformly for jc € and assume that is of 

ou 

constant sign for all (x, u) G £2 x (K \ {0} ) . Let /, g G W ~ 1 ' 5& (£2 ), f , g' G W 1 ~ p ' p id £2 ) with / > g 
in Q. (in the sense of distributions), /' > g' on d£2. Let «,v£ W l ' p (Q) solve (in the weak sense) 

-A p u = P(x,u)+f(x), -A p v = P(x,v)+g(x) in£2, 
u = j , v = g on a LI . 

Then one is interested in the following comparison results: 
(WCP) Weak Comparison Principle: Is it true that u > v in £21 

(SCP) Strong Comparison Principle: If u,v 6 c & 1 (£2), u ^ v, u > v in £2, is it true that u > v in £2 and 
7n ^ 3n ( x °) f° r an y x ° e ' Here, n is the unit outward normal to i2 on <9f2 . 

The Weak Comparison Principle (WCP) holds when < for which we refer to Tolksdorff I12T1 . 

ou 

dB 

The Weak Comparison Principle also holds when ^— > under the following assumptions and for 

ou 

Dirichlet boundary data: 

(A-l) > V (x,u) e £2 x (R\ {0}), j3(jc,0) > Vx G i2. 
(A-2) The problem 

— A p u = B(x,u) + f in £2 , 
m = on <5i2 . 

(where / G L°°(£2), f > in £2) admits a unique non-negative solution u G W l p (f2). 
(A-3) f,g^U°(£2),0<g<fon£2 and = = /' on d£2. 

This result is proved in @. However, for our purposes the zero Dirichlet data assumption in 
(A-2) and (A-3) is too restrictive. We show that this result also holds for inhomogeneous Dirichlet 
boundary data, that is, by relaxing the condition (A-2) and (A-3) to (A-2') and (A-3') respectively: 

(A-2') The problem 

—ApU=B(x,u)+f in £2, 
u = f on d£2 . 

(where / G L°°(£2), f > in £2 and f > on d£2) admits a unique non-negative solution u G 
W l -P(£2). 

(A-3') f,geL°°(£2),0<g<f on £2 and < g' < f on d£2. 



We prove the following results along the same lines as in (6). 
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Theorem 3.1 Let the assumptions (A-l) (A-2') and (A-3') hold then the WCP holds for bounded 
solutions. 

Proof. Let us denote = {h € L°°(i2) | h > in Q}. Given / e L+(Q) and /' e wW'^dfl) 

with /' > on dQ, define the nonlinear operator Tjp on L™(X2) by letting Tf p(u) = v, where v is 
the weak solution of 

-A p v = P(x,u)+f inn, 
v = f on d£2 

Since > and u > 0, it follows that /3 (x, u) > /3 (jc, 0) > 0. So, the right hand side in (13.21 i is 
non-negative as also the boundary data. By appealing to the WCP proved by Tolksdorff Ell we 
conclude that indeed Tt p(u) = v > and Tt f maps L™(£2) into itself. 
Claim. Let fi, fi,u\,U2 EL™(£2). If /i < /2, «i < «2 and /[ < /, then 7V /-'(hi) < T^fiuj) 

dp 
du 

Let vj = ^//(mi) and v 2 = T f2 ^(u 2 ). Then 



Indeed, following the condition > we conclude that /j* := j8(jc,Mi)+/i < j5(x,U2) + f 2 ='■ /| ■ 



-4pVi = /*, -A p v 2 = fl in f2, 

n=fi, V2=fi on dQ.. 

So, again by the weak comparison result proved in lED we obtain Vi < v 2 in £2. This proves the 
claim. 

Now, let u,v be bounded solutions of the non-linear pdes in d3.lt . To begin with, Tt f(u) — u and 
T g g i (v) = v. Now, using the claim we obtain the inequalities, 

< T fJI (0) < T fJ , («) = «, < T ?y (0) < T g j (v) = v 

We can then show by an inductive application of the claim that following chains of inequalities hold 

< T fif , (0) < Tjj (0) < • • • < T n S fl (0) < • • • < u = T fJ , (u) (3.3) 
< T g j (0) < Tl gl (0) < • • • < T» g , (0) < • • • < v = T gig , (v) (3.4) 

The pointwise limits m*(x) = lim 7? f ,(0) (*) and v*(x) = lim 7 1 " ,(0) (x) exist and must 

clearly satisfy TV /'(«*) = m* and T ? y (v*) = v* respectively. So, by the uniqueness assumption in 
(A-2'), it follows that u* = u and v* = v. 

Again, by applying the claim above, for any n > 0, we obtain T n , (0) < 7? « (0). Therefore, upon 
taking the limit as « goes to infinity we obtain v < m. This proves the theorem. □ 
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4 Existence and uniqueness for a nonlinear Dirichlet problem 

Let Ai be the first eigenvalue of the Dirichlet p-Laplacian as in ( 11.4b on a bounded domain £2. Let G 
be an open proper subset of £2. We prove the existence and uniqueness result for a nonlinear partial 
differential equation on & given Dirichlet data f > on d£2 . This shall be needed for applying the 
comparison principle of the previous section, later in Section|6] 

Proposition 4.1 Given f 6 W l ~ i i>' p {d 6") and f > on dff, the problem 

-A p w = X l \w\ p - 2 w inff, 1 
W = f on 60. J 



(4.1) 



admits a unique non-negative solution. 
Proof. Let us first prove that if a solution exists then it is non-negative. Let u be a solution of 



the above problem. As u > on dff, we obtain that u G Wq' p (&). Therefore, taking u as a test 



From this we obtain 



function, we have 

\¥u\ p ~ 2 (Vm, Vm~) dx — X\ \ \u\ p ~ 2 uu~dx. 

Jff 

[ \Vu-\ p dx = Xi [ \u-\ p dx. 

Jff Jff 

We cannot have u~ ^ 0, for otherwise, from the variational characterization of the first eigenvalue 
we can conclude that X\{ff) < X\ = Xi(Q). However, this cannot happen, being a proper open 
subset of £2 we must have X\(Q) < Xi(&). 

Existence. We denote by /' again aW l ' p (0) function whose trace on d & is /'. We can then obtain 
a weak solution of ( 14. \\ by minimizing the functional J(w) = Jg- Vw\ p dx — X\ (£2) j ff \w\ p dx on the 
affine space A ;=W^ p {ff)+f. Indeed, if w is a minimizer of J then we shall have 



dt 



J(w + t<j))= f |Vw|^ 2 (Vw,V<p)dx-Ai / \w\ p - 2 w<pdx Vq> (4.2) 

; = Jff Jff 



which is just the weak formulation of ( 14.1b . As A is a closed convex subset of the reflexive Banach 
space W l,p (iff), for showing the existence of a minimizer of J on A, it is enough to prove that J is 
coercive and weakly sequentially lower semi-continuous on A. 

/ is weakly sequentially lower semi-continuous on A: This is true since J n \ Vw\ p dx is lower semi- 
continuous for the weak topology on W ,p {ff) and J n \w\ p dx is continuous for the weak topology 
on W Lp (ff) due to the compact inclusion of W l ' p {ff) in L p {6). 

J is coercive on A: Let w n :— f + <p„ £ A be a sequence such that ||w„|| w i.p(^) — > °° as n — >• «>. If 
Iff \ w n \ p dx is a bounded sequence, then the coercivity is immediate. 



8 



Anisa Chorwadwala, Rajesh Mahadevan 



So, let us assume that j ff \ w n \ p dx — > oo as n — > oo. We may write w„ := f + <p„ with <p £ Wq ,p (&). 

Let B n := ^ ^"J . It can be argued, using the triangle inequality, that { \<p n \ p dx — > °° and 

ieWn\ p dx 

B„ — > 1 as n — >• oo. From the Poincare inequality on ^, we conclude that | Vq> n \ p dx — > oo as n — >• oo. 

Setting A„ := ^ " , we obtain using the triangle inequality, that jff\Vw n \ p dx —> °° and 

Jff v* < Pn\ p dx 
A n — } 1 as n — > oo. Now, 



where the last inequality has been obtained by applying Poincare inequality in the domain 6 . 
Since we have < k\(£2) < X\{0), since A n and B„ converge to 1 as n — s- oo, it follows that 
/ B n k\(Q.)\ 

A n [ 1 — — is bounded below by a positive constant C > 0. Once again, we have the coer- 



civity of /. 

Uniqueness. Suppose u,v are two different solutions of (14. 2\ in A. Let w\ :— Vlogw and w 2 '■ = 
Vlog v. As f(x) = \x\ p is a strictly convex function we have 

\w x \ p > \w 2 \ p +p\w 2 \ p ~ 2 (w 2 ,w 2 ~w l ) (4.4) 

and equality holds if and only if w\ — w 2 . If we prove that w\ — w 2 then we are done because in that 
case we will have = Vlogw — Vlogv = Vlog (-). That is, log f-J = k for some constant k. As a 
result we get u = e k v. But as u = v = f ^0 ond& we get u = v in & . Therefore, it suffices to prove 
that 

\w\\ p = \w 2 \ p + p \w 2 \ p ~ 2 {w 2l w 2 - wi) . (4.5) 

The proof of ( 14.5b is the same as the proof of Lemma 3 . 1 in Lindqvist [ 1 5 1 . We include the proof here 
for completeness. The function u solves ( 14.11) . We use u — v p u l ~ p as a test function in the equation 
for u. Similarly, we use v — u p v l ~ p as a test function in (14. U with v as a solution. Then we integrate 
by parts and sum the two identities. This new identity can be reduced to 

0= / Iwi \ p — \w 2 \ p — p \w 2 \ p ~ 2 (w 2 , w 2 — Wl)} 

Jff (4.6) 

+ v p {\w 2 \ p - |wi | p - p |wi \ p ~ z (w 1 , wi - w 2 )}\ dx. 

by using the following: 

V(u - v p u l - p ) = {l + (p - 1) (-)'} Vu - p (£j P ^ Vv, 



and, 



v(v-«V-O = {i + 0'-i)(") / '}vv-p(^ v« 
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But by ( 14.41 ) the integrand in ( 14.6b is non-negative (being the sum of two non-negative terms) and 
so, it follows from (14.6b that this integrand is equal to zero almost evrywhere in G. Therefore, each 
of the terms in the integrand must be zero. This proves (14. 5) . □ 



5 Shape derivatives of torsional rigidity and eigenvalue functionals 

Let Q be a smooth domain in M. N and let @ be a domain such that Q t C @, for t sufficiently small, for 
the smooth perturbations <P, associated to a smooth vector field V. Consider the Dirichlet boundary 
value problem on £2 t : 

—A p u =1 in Q t , | 
u = on dQ t . J 



(5.1) 



Lety,e < af 1 > a (A)bethe unique solution of problem (15.1b . Throughout this section y — y(Q) denotes 
the unique solution of (15. lb for f = 0. Denote (y r o 0,) \ a by / (f S R). We also denote the torsional 
rigidity E(A) by£(f). 



Proposition 5.1 The shape derivative of the torsional rigidity functional E(£2 t ) exists at t — and 

E{Q t ) = 



d_ 

dt 



dn 



dy 
dn 



p 

{V,n)dS. (5.2) 



r=0 

(Here, n denotes unit outward normal on dQ.) 

Proof. Let y be the unique solution of ( 15. Il l on £2 corresponding to t = 0. 

Step 1 : We first show that y' — > y strongly in W Q 1 ' P (X2). 

This can be obtained using the F-convergence (cf. Attouch [3 1, Braides [5 1, Dal Maso [7 1) of a 
suitable family of functionals. Consider the following family of functionals defined over Wq' p (£2): 

F(t,y):=- f \B;(Vy)\?Ytdx~ f yy t dx (5.3) 
p Jn Jn 

Since B* converges uniformly to / and y, converges uniformly to the constant 1, it is classical to show 
the F-convergence of the family of convex integral functionals F(t, •), as t — > 0, to the following 
functional 

1 



F(y):=F(0,y) = - \Vy\ p dx- y dx (5.4) 
p Jn Jn 

See Theorem 5.14 in Dal Maso |7| for instance. Furthermore, the family F(t,-) is equicoercive 
following the inequalities (12.7b and (12.5b . By standard results on F-convergence (cf. Theorems 7.8 
and 7.12 Dal Maso [7]), the minimizer of F(t , •) converges weakly in W ]' P (Q) to the minimizer of 
F(-) and the minima converge. Now, for each t e R, y t satisfies the equation 

[ \Vy t \P- 2 (Vy,y W )dx= f ydx V y/ e <C (A ) . (5.5) 
Jn, Jn, 
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By the change of variable <P, : Q — > Q t , the equation ( 15.51 ) can be re-written as 

|Z?;(V/)r 2 (A r (V/) ) V( P )dx = Jy t <l>dxV( P £ (5.6) 
Therefore, y' satisfies : 

-div(|B;(Vy)|"- 2 A f (V/)) -y t = in Q, } 



y'=0 ondQ. 



(5.7) 



which is the Euler equation for the minimization of the convex functional F(t,-) and therefore, / is 
the minimizer of F(t , •). Whereas, y, being the solution of problem ( 15. U for t = 0, is the minimizer 
of F . So, by the F-convergence result, we have the convergence of the minimum values 

lim^ o f--0 / \B;(Vy')\ p Y l dx= (--l) f \Vy\ p dx. (5.8) 



P J J a 

and the weak convergence in W ' ' p (Q ) , as t — > of to y. It remains to show the strong convergence. 

Since B* and y, converge uniformly to / and 1 respectively, and y t remains bounded in W ,P {Q), 
we can conclude from (15.8b that 



lim^o / \Vy , \ p dx= [ \Vy\ p dx. (5.9) 
Jn Jn 

Therefore, since the LP norm is uniformly convex, we can conclude from the weak convergence of 
Vy to Vy in LP (Q ) and the convergence of their norms (15. 91 that the convergence of V/ to Vy is 
strong in L P (Q). By Poincare inequality, as the LP norm of the gradients is an equivalent norm on 
Wq' p (Q), we obtain the desired conclusion. 

Step 2 : We observe that the torsional rigidity E(t) of the domain Q, is given by 
E(t) = -?— sup \ f $ 7t dx-- ( \B*(Vy)\Py t dx\ = ^— sup (-F(tJ)) (5.10) 

and the supremum is attained at <p — y' for / = y t o 0(f) and y, is the solution of ( 15.11 ) on i2,. 

Indeed, the supremum in the above corresponds to the negative of the infimum in the following 

inf (-/ \B*(W(p)\ p y,dx- [ <py t dx 



inf < — f \V(p\ p dx— [ (p dx\ 
W^iiiA I P Jn, Jn, J 



<t>eWg' p (n t ) <P Jn ' 

and this is attained by y t which is the solution of the Euler-Lagrange equation (15.11 ) on £2,. We can 
calculate this value which turns out to be 

F(t,y<) = - P — [ y t dx = -^-E{t) . 
P Jn, p 
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This proves our affirmation. 

E(t) -E(0) 

Step 3 : We now show that the shape derivative exists, that is the limit, lim,^o , exists 

and 

E(t)-E(0) p dF , . 

lim^o V ' y 1 = ^^- 0,y (5.11) 

t p — 1 at 



We obtain from the variational characterization (15.10) of E(t ) that 

E(t)-E(0) > -^— x (F(0,y)-F(t,y)) . (5.12) 

Thus, 

IkM^M > ^ Hm W^M = _^_^( 0>y ) . (5.13) 
40 t ~ p-l 40 f p-l dt 

Once again by applying the variational characterization of £(?) we have 

p — 1 

Therefore, by applying the integral form of the mean value theorem in the above in the first variable 

E(t)-E(0) < p F(oy)-F(*,y) 



p-l ? 

p /- 1 dF(st,y t ) 



p-l Jo 

In order to conclude the reverse inequality 



tasup £fc£ffi<_^^ (0 ,,) (5,4) 
40 f p-l dt 



it is enough to show that 



dF dF 

liminf — (st,y') = — (0,3/) for every s e [0,1]. (5.15) 

By a straightforward computation it is seen that, for any y S W 1 ' P (X2), we have 

^(s, W ) = i-(lf \B*(V V )\Py s dx- f yjsdx 
as os \p Jn Jn 

( -y s \By ¥ \Pdx+ [ \Vy\P- 2 ((B*)'Vy,ByY)dx- ( y> s ydx 
Jn p Jn Jn 

So, in particular, by taking y — y' we get 

?£(st,y t )=U~Y st KVyrdx+ J n \Vy'r 2 {(B:,yVy',B: t Vy')dx~ Jj.yd^j . 



12 



Anisa Chorwadwala, Rajesh Mahadevan 



Due to the strong convergence of V/ to Vy in L P (Q) and the analyticity of B s , y s in s, it is now 
straightforward to pass to the limit as t —> and we obtain easily, using (12.3b and (12.4b . that 

liminf -^(s*,/) = f / -divV |Vy| p c/x+ / |Vy| p - 2 (-(DV)*V3/, VyW* - / divVydx 
40 as \Jq p Jq Jn 

- £<M 

for every 5 G [0, 1], proving the claim (15.15b . 

Step 4 : To obtain the expression for the shape derivative (15. 2b . it is enough to integrate by parts in 

d F i 1 

the term — J n divVy dx which appears in the expression for -^-(0,y) = — — ^'(0) and apply the 

Pucci-Serrin identity ( 12. 9b . □ 



We now recall the shape derivative for the eigenvalue functional. Consider the eigenvalue problem: 

— A„u — \\u\ p ~ 2 u in Qti ) 

n (5 - 16) 
u = on dL2 t . I 

The first eigenvalue X\(t) :=X\ (Q t ) is simple and is characterized as the minimum of the problem 



Ai(A) :=inf< 



<pew { l' p (n t )\ . 



(5.17) 



We fix y\j .=yi (S2t) to be a corresponding eigenfunction which is positive (using the Rrein-Rutman 
theorem) and normalize it to satisfy 



/ \yu\ p 



dx — 1 . 



(5.18) 



For t = 0, we denote the corresponding eigenvalue and eigenfunction by Ai and y i respectively. 
Proposition 5.2 The map 1 1 — > X\ (t ) is differentiable at and 



A((0) = -(p-l) 



/ 


dy\ 


Ida 


dn 



(V,n)dS 



(5.19) 



Proof. As we have mentioned before, this result has been shown previously by de Lis and Garcfa- 
Melian [10] and by Lamberti [14] (also see Ly lfl6l ). This can also be proved along the same lines 
as in Proposition |5T| □ 
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6 Main Results 

Let < ro < r\, B\ be the ball B(Q,r\) and let Bo be any open ball of radius ro such that Bq C B\. 
Consider the family & = \B\ \Bo\ of domains in M. N . We study the extrema of the functionals 
E(£2) and X\ (£2) over associated to the problems ( II. U . ( 11.3b respectively. 
We state our main results : 



PutA) = fi(0,n)\fi(0,r ). 

Theorem 6.1 77ze minimum value of the torsional rigidity functional E(£2) on & is attained only 
when £2 — £2q> i.e., when the balls are concentric. 

Theorem 6.2 The first Dirichlet eigenvalue X\ (£2 ) is maximum on & only when £2 = Qq, i.e., when 
the balls are concentric. 

Before proceeding to the proof we make the following observation and reduction. The functionals 
to be optimized are invariant under the isometries of M. N . Therefore, it is enough to study these 
optimization problems for the class of domains £2(s) :— B\ \B(se\,ro), < s < r\ — ro where e\ 
is the unit vector in the direction of the first coordinate axis. In order to study the optimality of 
the domain £2{s) in the class & we need to study perturbations of the domain which correspond 
to translations of the inner ball along the direction of the first coordinate axis. For this purpose we 
consider a smooth vector field V(x) = p(x)e\ VxSBi where p : R N — > [0, 1] is a smooth function 
with compact support in B\ such that p = 1 on a neighborhood of B(se\,ro). Let {<£(}, g r be the 
one-parameter family of diffeomorphisms of B\ associated with V. We see that, for t sufficiently 
close to 0, <P t (Q(s))) = Q(s + t ). So, if we define : (ro — r\,r\ — ro) — > M as follows: 

j(s)=E(n(s)) and j^s) = Ai(fl(*)) 

we see that the minimization of E in the class & corresponds to studying the minimum of j on the 
interval (ro —r\,r\ — ro) and that the problem of maximization of in the class & corresponds to 
studying the maximum of j\ on the interval (ro — r\,r\ — ro). Also, the shape derivative of E and A[ 
at Q (s) for the vector field V are the ordinary derivatives at s of j and j\ respectively. We have seen 
in Proposition 15 . 1 1 and 15 .21 that these shape derivatives exist and so the derivative of both j and j\ 
exist. The optimization problems can be studied by analyzing the sign of the derivatives of j and j\ . 

First, we note that both j and j\ are even functions and since they are differentiable, we have 
/(0) = = / 1 (0). 

We shall adopt the following notations. Given s in (0,ri — ro) we simply denote Q (s) as Q and 
B(se\ , ro) as Bq and n shall denote the unit outward normal to £2 on d£2 . Let H denote the hyperplane 
H := {x = (x\ ,X2, ■ ■ ■ ,xn) £ M. N \x\ = s}. Let r# be reflection function about H. We define to be 
the subdomain G := {x £ Q \x\ > s} in Q. Then we see that the reflection of ff about H, namely 
&' = rn ( &) is contained in B \ , whereas Bo is symmetric with respect to H. Thus &' C Q . For x£ & , 
let x' denote the reflection of x about H, namely, the point r#(x). With these notations, if y be the 
solution of the equation ( 11.11 ) in £2 then, from the expression of the shape derivative (15.2) for E we 
obtain that 

dy 



xedB 



on 



p 

n\(x) dS (6.1) 



14 



Anisa Chorwadwala, Rajesh Mahadevan 



since V is zero on dBi and since, for all x 6 SBq, <V,n> (x) = pe\ -n(x) = n\ (x), the first compo- 
nent of the normal vector. Similarly, if yi be the solution of ( 11.3b in £2 then, from the expression of 
the shape derivative of X\, viz. ( 15.191 l, we obtain that 



j[( s ) = -(p-l) [ 



dn 



«i(x) dS 



(6.2) 



Proof of Theorem 16.11 Let y be the solution of the boundary value problem (11. Q in Q . We recall 
that y £ < rf? l ' a (Q) by regularity results in Tolksdorff 11221 . and by the strong maximum principle (cf. 
Theorem 5, Vazquez [23]) we have y > in Q. We now consider the subdomain 6 and let us define 
y on by y(x) :—y(x') the value of y at the reflection x' of x about H. Let us note that || (jc) = ^ (x') 
and «i (x') = — «i (x) for all x £ 3Bq. Now, we may rewrite the expression (16.lt as follows: 



f(s) 



■edB ndff 



edB nd& 



on 
an 



n i (x) rf5 - 



an 



x'edB ndff' 



n\(x) dS. 



an 



n\(x ) ofS 



(6.3) 



We shall show that /(s) > V 5 € [0,ri — ro) and is zero only if s = 0. We have already observed 
that /(0) = by symmetry considerations. It is clear that «i (x) < Vx e d& n dBo<~)H c . So when 
s ^ 0, we shall prove that / (i) > by showing that 



dy 
dn 



on 



(x)<^(x)<0 Vx<Ed<ffr\dB nH c . 



(6.4) 



We shall prove inequality ( 16.4b in a few steps. 

Step 1 : First we prove that — < on 3Bq. 

on 

We begin by noticing that at every point xq on <9Z?o, the interior sphere property holds, that is, there 
exists an open ball B = Br(zq) C £2 such that dB n 3Bq = {xq} and the unit outward normal n 
to Q and to B coincide at xq. For K > and a > we define a function b : B — > K as b(x) = 
^ (e-"l*-*ol 2 _ e -«« 2 ^) . We have b > in B r (z q )\Br (z ), b(x ) = 0, in fact, ft = on dB and that 

— (xo) < 0. Moreover, it can be shown that 
dn 

-A p b(x) - -2Ke- a \ x - z ^\Vb\ p - 2 (2a 2 \x-z \ 2 ~Na) . 
We may therefore, choose a large enough (independent of K) so that 

-A p b <0 in B r {zq)\Br(zq). 



(6.5) 



We know that y satisfies (II . lb . Since y is bounded below by a positive constant on dBn (zq), we may 
choose K small enough so that b < y on dBit (zo). Thus we have 

b<yond(B R (zo)\BR{zo)), -A p b<0 and - A p y > in B R (z ) \Br (zq). (6.6) 
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Then by the WCP of Tolksdorff fl2T| we have b < y in B R (z ) \Br (z ). This, along with b{x ) = = 

db . dy , db dy . 

y(xo), implies that — (xq) > — (xo). Since — (xq) < we get — (xq) < 0. We have thus proved 
on on on on 

dy 

that <0on dB o nd0T\H c . 
on 

Step 2 : Now we prove the first inequality in (16.4b . 
On ff, the function y satisfies 

-A p y =1 in &, 

y = ond&ndB , 
y = y * ond&nH, 
y = ond&ndBi; 

while y satisfies 

-A p y =1 in ^, 

y = ondffndBo, 
y = y* ond&HH, 
y>0 ond&ndBu 

where y* denote the common value of y and y on H. Therefore, we have 

-A p y = ~A p y in ff, 
y>y ondff. 

So, by the WCP of Tolksdorff [21] we get, y > y in & . Since, y = = y on 3Bq n dff we conclude 
that 

-/W<^Wforall^e5B n5^. (6.7) 
on on 

By the result of Step 1 and (16.7b . we can obtain a neighborhood jV of <9Z?o H <9 and positive 
numbers 77 , £0 small enough so that 

l*Vy(x) + (l -t)Vy(x)\ > 77, ty(x) + (l-t)y(x)<£ Vt E [0, 1], M x £ ~JY . (6.8) 

Let w:=y-y, then w>Q on dff with w = on d@f\ (HUdBo) and w > on d&ndBi. Let 
A : R N -t R N be the map as defined in (l210t . Then, we have 

div(A(Vy) -A(Vy)) = 4 /3 y - A p y = Qm,yV. 

By Mean Value Theorem we get, 

tj^.( / 1 (VA ! )(fVy+ (1 -t)Vy)dt,Vw) = 0. 
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Let oy(jc) = Jo ^(tVy(x) + {l-t)Vy(x))dt. Then, w satisfies 



* d_ ( dw_ 



m.vV, 



w>0 in^K, 

w = ond,yV ndB Q . 



(6.9) 



By dZTTT ). [fly(x)]f. =1 > min{l,/7 - 1} /J |fVy(x) + (1 -f)Vy(x))|'>- 2 I if and so, by (EHJ we con- 
clude that [fly (*)],-/_! is a uniformly positive definite matrix when * S =yT. Therefore, by the maxi- 
mum principle for uniformly elliptic operators (cf. Theorem 5, Ch. 2, Protter and Weinberger 11171 ). 
since w is a non-constant function, it follows that the minimum of w will be attained on djV . Since 
wii^jv w = it follows that w > in Further, by the same argument as in the Hopf 's Lemma for 

uniformly elliptic operators (cf. Theorem 7, Ch. 2, Protter and Weinberger IfTTl ) we have -z— < 

on 

on djV n <9Bo- That is, V.iG d0ndB o we have the following: 

Iw<|w<o- 

□ 



Proof of Theorem 16.21 Recall from Section [5] that y\ is the principal eigenfunction of (11.31 > (that 
is, the unique solution of ( 11.3b for A = Ai (i2 )) characterized by y i > in i2 and y p x dx~ 1 . We 
now consider the subdomain and let us define y\ on <?? by y\{x) \— yi(x r ) the value of y\ at the 
reflection x' of x about //. Let us note that ^{x) = ^{ x ') an d n i( x ') = — ii(*) for all jc 6 <9Z?o- 
Now, we may rewrite the expression (16.21 as follows: 



/ { 


dyi 


(x) 


JxedB ndff i 


dn 


f I 


dyi 


(x) 


JxedB ndff ( 


dn 



p , 

n\ (x) dS — (p — 1) / 

Jx' 



an 



edB ndff' 
n\ (x) dS. 



on 



ni(x') dS 
(6.10) 



We shall show that < V s £ [0,n — ro) and is zero only if s = 0. We have already observed 
that /j (0) = by symmetry considerations. It is clear that n\{x) < V x G d & n 3Bq n H c . So when 
s ^ 0, we shall prove that /j (s) < by showing that 



on on 



(6.11) 



We shall prove inequality ( 16.1 U in a few steps. 

Step 1 : First we prove that < on <3Bo- 
on 

As in the proof of Theorem l6.1l we begin by noticing that at every point xq on 3Bq, the interior sphere 
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property holds, that is, there exists an open ball B = Br(zq) C Q such that dB D 9Bq = {xq} and the 
unit outward normal n to £2 and to B coincide at Xq. We recall that y\ satisfies (11.3b and y\ > in 
Q . We construct an auxiliary function b, as in the proof of Theorem 16.11 with a sufficiently large 

so that — A p b < on Br(zo) \ Br (zq) and K small so that b < y\ in d (Br(zq) \Br(zq)^, to obtain 

9y\ 

b < yi in Br(zo) \Br (zq) and consequently to obtain < on 9Bq. 

7 . . an 

Step 2 : Now we prove the first inequality in (16. lib . 

On G, the functional satisfies 



whereas y~\ satisfies 





= Xiy[ 


in 6, 




y\ 


= 


on 96 


'ndB Q 


y\ 


* 

=y\ 


on 96 


'f)H, 


y\ 


= 


on 96 


'ndBi 


A P h 




in ff, 




y\ 


= 


on 96 


'ndB 


h 


=y\ 


on 96 


'f)H, 


h 


>o 


on 96 


'ndBi 



where y* denotes the common value of y\ and y\ on H. Therefore, we have 

-A p yi = hff 1 in &, 
y\ > yi on 96". 

Therefore, by the WCP (cf. Theorem 13.11 proved in Section [3} for ^>0we get, y~\ > y\ in 6. 
Since, y\ = = y\ on 9Bq n 9& we conclude that 

^"W < ^-W forallxG 5B na^. (6.12) 
on on 

By the result of Step 1 and ( 16.12b . we can obtain a neighborhood JV of 9Bq n9G and positive 
numbers rj , £o small enough so that 

\tVyi(x) + (l -t)Vyi(x)\ > 77, tyi(x) + (1 -t)yi(x) < Bq Wt e [0,1], VxG^F. 

Let w :=yi -yi, then w > on with w = on dt?n(HUdB ) and w > on 90 ndBi. We 
have 

-</fv(A(Vyi) -A(Vyi)) = -4 p yi +^1 - h (yV ~yV) > Oin ./K, 
where the map A is as defined in (12.10b . By Mean Value Theorem we get, 

-£^^V^)(»Vyi + (l-0Vyi)*,Vw^ >0. 
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Let fly(jc) = Jq gjf(fV;yi(jc) + (1 -f)V;yi(jc))df, then w satisfies 



(6.13) 

w>0 in^K, 



vr 



= ond<j¥ ndB Q . 



As in the proof of Proposition 16.11 we observe that the matrix [fiy(x)]yj =J is a uniformly positive 

definite matrix when x G JV . Then by the maximum principle for uniformly elliptic operators (cf. 
Theorem 5, Ch. 2, Protter and Weinberger flTT]), since w is a non-constant function, it follows that 
the minimum of w will be attained on d^V. Since infg ^ w = it follows that w > in ^/K. Further, 
by the same argument as in the Hopf's Lemma for uniformly elliptic operators (cf. Theorem 7, Ch. 

2, Protter and Weinberger [17]), we have < on dJV ("1 dB Q . That is, V x G 5 ^ n <?fio we have 
the following: 

< -5-w <o. 

OB (7« 



□ 
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